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Abstract. We discuss the quasi-adiabatic approximations to the three-body 
wavefunction in breakup processes, clarifying the assumptions underlying the 
model. This suggests alternative approximation schemes. Using different the- 
oretical three-body models, calculated differential cross section angular distri- 
butions for the 11 Be(p,d) reaction, for which new preliminary data have been 
reported at 35 MeV, are presented. We show that calculations are sensitive to 
the inclusion of deuteron breakup and to the breakup model used, particularly 
if used to deduce absolute spectroscopic information on the + and 2 + 10 Be 
core state parentages. There is also considerable sensitivity to the model used 
in calculations of the relative cross sections to the two states. 

1 Introduction 

In nuclei near the dripline, the separation energy of the last nucleon(s) becomes 
extremely small. Compared with the common 6-8 MeV in stable nuclei, many 
dripline nuclei have either a nucleon or a two-nucleon separation energy that 
is less than 1 MeV. The neutron-density distribution in such loosely bound 
nuclei shows an extremely long tail, called the neutron halo. Although the 
density distribution of a halo is very low, it strongly affects the reaction cross 
section and leads to new properties in such nuclei. The study of halo nuclei 
is interesting as they involve new structures and surface densities, and thus 
provide a stringent test of theoretical models of nuclear structure and reactions. 

Recently, the use of low energy single nucleon transfer reactions for structure 
studies of exotic nuclei have attracted attention fl]]-||. Because of the simplicity 
of the theoretical interpretation of these reactions, they are thought to provide 
an important source of the information about the structure of halo nuclei. 

Single nucleon transfer reactions, such as the (d,p) and (p,d) reactions, have 
been a reliable tool in nuclear spectroscopic studies of stable nuclei, determin- 



* Web address: 



bttp: / /www.gantcp.cdu.tr^ayilmaz 



2 



ing positions, spins and parities of nuclear states. Developments in radioactive 
beams technologies are now producing beams of nuclei at and near the neutron 
and proton driplines, including neutron halo states with very weak binding. 
The positions and the ordering of the nuclear single particle levels in such sys- 
tems, with large neutron excesses, has yet to be clarified. These exotic systems 
typically have no bound excited states and so traditional spectroscopic meth- 
ods are inapplicable. Using the single nucleon transfer process, on a deuteron 
or proton target, in an inverse kinematics experiment, is then an attractive - 
although still very difficult - alternative tool for mapping out nuclear structures 
near the driplines. 

11 Be is a good example of such a system having a single neutron separation 
energy of only 0.5 MeV. It is now understood that the 2si/2 neutron single 
particle state in this region is lowered and that a dominant component of the 
11 Be ground state is produced by the coupling of a 2si/2 neutron to a 10 Be (g.s., 
+ ) core; with a smaller but significant component in which a ld 5 / 2 neutron is 
coupled to a 10 Be (excited core, 2 + ). This nucleus is of particular theoretical 
interest because the ground state parity is exactly opposite to what one naively 
would expect from the spherical shell model. There have been many theoretical 
attempts to address the problem of parity inversion of the ground state and 
the first excited state of the 11 Be and evaluating the extension of its neutron 
halo. Most of them Q] correctly reproduce the parity inversion but make very 
different predictions about the degree of coupling of the 11 Be ground state 
with the first 2 + excited state of 10 Be at 3.368 MeV, the ratio of spectroscopic 
factors in these different models S(2 + ) / S(Q + ) varying from 0.07 to 0.73. In 
this context, the recent experiment at GANIL |l|] is expected to measure the 
transfer reaction cross section angular distributions at forward angles to the 
10 Be ground and 2 + states, with the expectation of clarifying quantitatively 
this 11 Be (g.s.) admixture, and the experimental determination of this ratio of 
spectroscopic factors by means of a neutron pick-up reactions. 

The purpose of the present work is to investigate the importance upon 
such transfer reaction spectroscopic studies of the inclusion of the deuteron 
breakup degrees of freedom in the theory used to analyze measured cross section 
observables. Of importance in analyzing the experimental measurements will 
be the extent to which the magnitudes of the calculated cross sections, and 
particularly the ratio of the cross sections to the ground state and 2 + core 
states of 10 Be, are affected by the inclusion of three-body channels. 

We discuss the 11 Be (p,d) reaction within different theoretical three-body 
(n+p+ 10 Be) models. We calculate the transfer amplitude using the prior form 
of the (p,d) matrix element, thus the transition interaction is the n-p interaction 
and we need a full (three-body) description of the n+p+ 10 Be system in the final 
state. For the description of this final state we have used 

i) the adiabatic (AD) model ||, 

ii) the quasi-adiabatic (QAD) model and 

iii) the projection operator (POA) approach ||. 
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To clarify the importance of the breakup corrections to the transfer cross 
sections, we also perform Distorted- Waves Born Approximation (DWBA) cal- 
culations. For the ground state transition a preliminary test of the models 
against more exact Coupled Discritized Continuum Channels (CDCC) model 
calculations has been carried to provide an assessment of the reliability of 
the three-body models used at the present incident energy. 

In Section 2, we also present a re- formulation of the quasi-adiabatic model 
for nuclear reactions which makes clear the approximations inherent in the 
model. We show that the so called quasi-adiabatic breakup wavefunction also 
contains the non-adiabatic corrections to the elastic part of the wavefunction. 
This was a significant uncertainty in the formulations given in Ref. §,@. In 
addition, in the same section, an alternative quasi-adiabatic scheme, the pro- 
jection operator approach (POA), based on the use of projection operators is 
introduced for a better description of the non-adiabatic elastic channel contri- 
butions. 

Current experimental activity in the area of light- neutron rich and drip- 
line nuclei now dictates the rapid development of calculable theoretical models 
for reactions and scattering of effective few-body systems. Though, the CDCC 
model has been spectacularly successful in understanding a wide range of data 
and phenomena in light- and heavy-ion three-body systems, and can provide 
benchmark calculations against approximate models such as AD, QAD and 
POA, it is unlikely ever to find an application to the solution of many-body 
problems with more than three interacting particles. Convergence problems 
may also become more serious when very weakly bound halo nuclei are involved. 
To date such systems have been analyzed using few-body Glauber based models 
[ fio| whose first step is an adiabatic treatment of the internal degrees of freedom 
of the projectile. However, the application of the technique at lower projectile 
energies, or to the case of including Coulomb breakup, has shown that the model 
begins to break down due to known inadequacy of the adiabatic approximation, 
see e.g. [D,Jll|], at such energies. This failure requires the treatment of non- 
adiabatic effects in the model. Quasi-adiabatic ideas are therefore an obvious 
and necessary generalization of these and non-eikonal adiabatic models. Hence, 
we first clarify the formulation of quasi-adiabatic models within the three-body 
context in the following Section. The calculation methods and the results of 
our calculations are discussed in Section 4 and 5, respectively. 

2 The Three Body Models 

The n+p+nucleus Hamiltonian will be written 

H = H np + T R + V p (r p ) + V n (r n ) + V c {r p ) (1) 

where H np = (T r + V np ) is the n-p Hamiltonian, Tr the center of mass kinetic 
energy operator and V n and V p the neutron- and proton-target effective inter- 
actions, regarded as local optical model potentials. Vc is the Coulomb field and 
is assumed to act on the n-p center of mass. 
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The total wavefunction ^ in the exit channel of the pick-up reaction of 
interest satisfies 

[E - H np - T R - V(r, R)] &(r, R) = (2) 



V(r, R) = V p (r p ) + V n {r n ) + V C {R) (3) 

where r is the relative coordinate of the n-p pair and R is the center of mass 
coordinate. 

In the context of (p,d) reactions enters the transition amplitude 

T pd = (&(r,R)\ V np \x ( ^(r P )Mrn)) (4) 

with (j) n the neutron bound state and X incoming proton wave function. 
Throughout this paper we restrict the formalism to S-wave n-p relative mo- 
tion for simplicity. In zero-range approximation then it is the wavefunction at 
coincidence, ^(O, R) which is of importance. 

Theoretical methods based on the adiabatic treatment of one or more quan- 
tum mechanical degrees of freedom have played an important role in nuclear 
physics. More specifically, the Johnson-Soper application of adiabatic ideas |^| 
to nuclear breakup effects played a key role in the development of models of 
breakup processes in the three-body systems. However, in the AD treatment 
it is assumed that the excitation of the projectile is to states in the low en- 
ergy continuum; its treatment of possible high energy breakup contributions is 
therefore naturally suspect. 

An improved treatment of these higher energy breakup configurations is 
provided by the QAD method calculations which takes approximate account of 
modifications to the center of mass energy of the n-p pair in breakup configu- 
rations through the use of a mean breakup relative energy for the continuum 
states. It thus breaks the degeneracy with the elastic channel. The inclusion of 
these higher energy breakup components via the QAD model calculations Jt| 
for large transferred angular momentum (d,p) transfer reactions (at energies 
E/ A « 40 MeV) led to significant deviations and systematic improvements over 
the adiabatic model in the description of experimental data. 

Although the quasi-adiabatic calculations produced an improved description 
of the measured observables, the theoretical justifications of the assumptions 
made in the model have not yet been studied. The work described in this paper 
is thus concerned also with the clarification of these theoretical uncertainties. 
In addition, we re-formulate the quasi-adiabatic theory to give a more general 
formalism, approaching the three-body problem in a different way. This alter- 
native formulation provides a clear understanding of the assumptions made in 
the original Q quasi-adiabatic theory. 

However, the quasi-adiabatic model does not include back-coupling mod- 
ifications to the elastic component of the wavefunction. For the inclusion of 
these modifications, we also develop here an alternative approximation scheme 
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for the treatment of quantum mechanical three-body systems using the unified 
theory of Feshbach. 

2.1 Adibatic Approach 

A considerable simplification in the solution of three-body equation is achieved 
using the Johnson-Soper adiabatic approximation JjjJ. The approximation, in 
which the dynamics associated with the continuum coupling collapses to an 
effective two-body problem for the n-p center of mass motion, involves the 
replacement of H np in the exact three-body Schrodinger equation, Eq.(2), with 
—Ed, the deuteron binding energy. Therefore we have 

[E c , n . - T R - U(r,R)] W AD {r,R) = (5) 

where E c , m , = E + Ed is the energy of the outgoing deuteron in the center of 
mass frame, and U(r, R) is the angle average of V(r, R) with respect to r. 

The assumption made is that the dominant breakup configurations are 
states of low relative n-p energy Ek when compared E with, such that little 
error is made by assuming 

E-ek = E + e d (6) 

In spite of the efficiency and success of the adiabatic approach, experimental 
transfer reaction data p"2[ clarified that some physical contributions are miss- 
ing from the calculation of the reaction amplitude. These involve the transfer 
of particles in large n-p relative energy configurations. As the adiabatic ap- 
proximation is formulated under the assumption of low energy n-p breakup, 
its treatment of high energetic contributions is naturally suspect. An improved 
treatment of the higher energy breakup configurations, such as quasi-adiabatic 
calculations, is thus required. 

2.2 Quasi-adiabatic Model 

Under the restrictions to S-wave relative n-p configurations, one separates ^ 
into its elastic and breakup parts. Thus 

V BU (r, R) = &(r, R) - W EL (r, R) (7) 

and therefore 

[E - H np - T R - U(r,R)] V BU (r,R) = [U(r,R) - U opt (R)] V EL (r,R) (8) 

where U° pt \l/ EL = [E c , m , ~T R ]^ EL . At this stage, the quasi-adiabatic model 
assumes: 

I. that H np is replaced by an average energy depending at most parametri- 
cally upon r, and 
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that the elastic part of the wavefunction generated by the adiabatic model 
■q;AD,EL j g an accura ^ e representation of >P EL . Thus the quasi-adiabatic ap- 
proximation to *]s BU solves 

[E-e-Tr- U(r,R)] V QAD > BU (r,R) = [U(r,R) - U AD < opt (R)] V AD > EL (r,R) 

(9) 

where e is taken JtJ as the expectation value of H np in \p AD ' BU in each partial 
wave, and u AD ' opt V AD > EL = [E c . m . - T R ] $ AD < EL . 

Thus, the quasi-adiabatic approximation removes the degeneracy of the n-p 
center-of-mass energy in breakup configurations by introducing a positive mean 
energy for the continuum states. We note that the model provides, however, no 
prescription for modifications to the elastic component of the wavefunction. 

2.3 An Alternative Formulation Of The Quasi-adiabatic Model 

We present here an alternative development of a formal theory for the quasi- 
adiabatic method, which clarifies that such a model can be introduced by mak- 
ing only one single approximation, unlike the formulations of the model in Ref. 
H , . Another advantage of this formulation is that one sees how to treat cor- 
rections to both the clastic and breakup components of the wavefunction, and 
to derive an iterative scheme for such changes. This was a significant uncer- 
tainty in the original quasi-adiabatic formulation of Amakawa et al. 0, in which 
it is stated that the elastic wavefunction is assumed unchanged, regardless of 
changes made in the breakup piece of the wavefunction. 

Under the restriction to S-wave relative n-p configurations, a formal devel- 
opment of the quasi-adiabatic theory proceeds by decomposing the projectile- 
target three-body wavefunction into the adiabatic wavefunction plus the cor- 
rection term, i.e., 

&(r, R) = <P AD (r, R) + AW(r, R) (10) 

where AP accounts for non-adiabatic corrections to the breakup and elastic 
channels, and has only outgoing waves since <P AD already satisfies incoming 
wave boundary conditions. Upon substitution in the Schrodinger equation then, 

[E - H np -T R - U(r,R)] AV(r,R) = (H np + e d )* AD > BU (r, R) (11) 

where the source term has infinite range. To proceed we use the inhomogeneous 
equation for \p AD > BU i 

V AD ' BU (r,R) = [Ec m . -T R - C/(r,i?)] _1 [U(r,R) - U AD > opt (R)] W AD ^ EL (r,R) 

(12) 

from which 

AV(r, R) = [E- H np -T R - U(r, R)] ~* {H np + e d ) [E c . m . — Tr — U(r, R)f3) 
x [U(r,R) - U AD < opt (R)] V AD < EL (r,R) 
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It follows that 

AW(r,R) = {[E- H np -T R - U(r,R)] ^ - [E c . m . - T R - U(r,R)] "f^) 
x [U(r,R) - U AD - opt (R)] & AD < BL (r,R) 

and hence we can write /^(r, R) = A<Pi(r, R) — A^ir, R) with 

[E - H np - T R - U(r,R)} A^(r,R) = [U(r,R) - U AD <° pt (E)] * AD ^ EL (r,R) 

(15) 

[E c . m . -T R - U(r,R)] AV 2 (r,R) = [U(r,R) - U AD '° pt (R)] V AD ' EL (r,R) 

(16) 

The equation for A\Pi leads to the original quasi-adiabatic wavefunction. 
A\p2, precisely the adiabatic breakup wavefunction and has no overlap with 
the elastic channel. Thus A&i must include both elastic and breakup non- 
adiabatic corrections. The elastic piece can be extracted by projection. This 
makes clear that the assumption *P EL » \p AD ' EL of the original formulation 
is unnecessary. Therefore the reduction of the exact three-body equation to 
the quasi-adiabatic model requires only the replacement of H np by an average 
breakup energy e. The details of the similar resulting partial wave expansions 
and solution of the equations can be found in Sec. Ill of Ref. @. 

2.4 Projection-Operator Approach 

In the previous analyses, the non-adiabatic elastic corrections are clearly 
treated incorrectly. They solve equations with the wrong asymptotic energy. 
In case of deuteron scattering the non-adiabatic elastic corrections are small 
but in general this may not be the case. A more correct treatment requires use 
of projection operators on P and off Q the projectile ground state wavefunction, 

W = PW + Q<f = f EL + W BU (17) 

The usual reduction of the Schrodinger equation can be shown || to yield 
coupled equations 

[E - QH np Q -T R - U(r,R)] <P BU (r,R) = [U(r,R) - U° pt (R)] V EL (r,R) 

[E c . m . -T R - U dd (R)] V EL (r,R) = \Mr)) (Mr)\ U(r,R) \* BU (r,R)) (19) 

with Udd the Watanabe potential and 4> d the deuteron ground state wavefunc- 
tion. Again at this stage H np is to be replaced with an average excitation energy 
e, in each partial wave. The equations can be solved by iteration || starting 
from the adiabatic model estimates. 
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3 Experimental Aspect Of The Reaction 

In the past, transfer reactions with light and heavy ions have been a major 
source of spectroscopic information in stable nuclei, including spin and parity 
assignments to nuclear levels, measurements of occupation probabilities and 
wave functions in the ground and excited states of the daughter system. Pro- 
vided that beams of exotic nuclei with high enough luminosity will become 
available it is tempting to use transfer reactions in the same spirit also for 
spectroscopic studies in exotic nuclei. The important difference to former fixed 
target experiments is the low intensity of secondary beams and the use of in- 
verse kinematics. Under these conditions only transfer reactions with rather 
large cross sections can be investigated. 

Since exotic beams are obtained as secondary beams they depend on the 
production rates for unstable nuclei in primary reactions with a stable beam. As 
a consequence, their intensities are reduced by several orders of magnitude with 
respect to the intensities of the primary beam. As a rough rule, the production 
cross sections for these beams fall by approximately one order of magnitude 
for each mass unit going further away from stability g . Thus, one is restricted 
to reactions with larger cross sections. Targets should be chosen as thick as 
possible and detectors with large solid angles and high detection efficiency 
should be used. 

Obviously, the abundance of data to which one is accustomed from work 
with stable beams can certainly not be expected for radioactive beams. In re- 
actions with inverse kinematics the cross-sections are focussed more in forward 
direction in the laboratory system. Thus, the full yield of the nucleus to be 
investigated is collected in a rather small angular cone and total cross section 
measurements should already be feasible at much lower beam intensities. 

As a representative example, we consider here an inverse kinematics reac- 
tion with beam of neutron rich nucleus n Be. The n Be (p,d) 10 Be reaction, 
leading to bound states in 10 Be, has been studied for the first time 0,0, using 
a secondary 11 Be beam of 35 MeV/nucleon. Angular distributions up to about 
15 degrees (in cm.) were measured by detecting 10 Be in a spectrometer and 
coincident deuterons in a position sensitive silicon detector array. Their prelim- 
inary analysis provides evidence for a large core excitation component in the 
structure of 10 Be ground state. 

4 Calculations 

In the present work, the 11 Be system is interpreted as a 10 Be core, with one 
excitable state, and a loosely bound neutron. We have only included the + 
ground state and 2 + excited state of the core. That means, in the ground state 
of 11 Be, the neutron oscillates between a 2si/2 and lds/2 orbital, exciting the 
core from a + to a 2 + state (and de-exciting it). 

Since a proper model provides modifications to the three-body deuteron- 
channel wavefunction, one expects the calculations resulting from this method 
to improve upon the predictions on reaction observables. Thus, in order to 
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establish the significance of the modifications introduced by the models, and 
for future comparison to the experimental data, we have evaluated differential 
cross-section angular distributions at 35 MeV incident energy for the 11 Be 
(p,d) 10 Be (g.s and 2 + ) transitions and the trends of these observables have 
been compared with each other. 

Three-body wavefunctions in the exit channel of the reaction are calculated 
using the adiabatic, quasi-adiabatic, projection operator and CDCC methods 
(the later using the code FRESCO |]|). The main objective of the CDCC 
calculations, which do not include spin-orbit interactions due to the difficulty 
in including this term in FRESCO, is to provide a critique of the theories 
implemented here, rather than to seek a realistic comparison with experiment. 
The preliminary test (for 2sj /2 transition) of the AD and extended AD models 
using the CDCC calculations, discretizing the n-p continuum into 6 bins from 
to 25 MeV excitation energy, provide an insight into their reliability at the 
lower incident proton energy of interest. 

The quasi-adiabatic and projection operator model calculations are iter- 
ated in the sense that the mean breakup energy for the continuum, £jl, in 
each partial wave is calculated from the latest best estimate of the breakup 
wavefunction, i.e. 



-(*) m _ ( il) *?F^R)\ H np |Wflffi (r,fl))„ 



where bra-ket denotes the radial integration over r and (i) represents the itera- 
tion number. Although this prescription reproduces meaningful breakup ener- 
gies at medium incident projectile energies for heavy targets ||, it breaks down 
at lower incident energies and cannot be iterated due to the calculated unphys- 
ical ejl values which are larger than the center of mass energy. To overcome 
this problem for the underlying reaction, we consider another formulation for 
the mean breakup energies. 

Starting with an exact definition of the continuum channel breakup wave- 
function in partial wave form , 

/■oo 

<P?U(r,R)= dkMr)XJLk(R) (21) 
Jo 

where </>fc(r) is a triplet n-p scattering state with asymptotic normalization 

(22) 

such that / °° drr 2 <j)k(r)<j)k' (r) — 5{k — k'), one can rearrange the mean energy 




prescription, Eq. (20), in the form 



sjl(R) 



(*?L J \H np \*fF) r 



JL 



U\ q,BU\ 



(23) 
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JT dk X * JLk {R)ekXJLk{R) 
Jo°° dk X * JLk (R) XJLk (R) 



r 



(24) 



where * denotes complex conjugate and H np (j)k(r) = £k4>k(r) in which e k — 

»- 2 .2 

One needs at this stage to evaluate the XJLk, by integration 



where &jl is calculated approximately by the QAD and POA theories. To be 
in consistent with the CDCC calculations carried out in this work, we set the 
maximum value of k for the integral in Eq. (23) to 0.78 fin -1 that corresponds 
to 25 MeV for the relative breakup energies. 

The calculated ejl values show significant changes from the zeroth order 
estimate based on \]/ AD ' BU but beyond the next iteration such changes are not 
reflected in changes in the calculated wavefunctions or in reaction observables. 

The calculated three-body wavefunctions using the same inputs, in partial 
waveform, provided by the present models are employed in a modified version 
of the program TWOFNR [jlij for the evaluation for the reaction observables, 
performing the zero-range approximation. The zero-range approximation is ex- 
pected to produce a more accurate model for in particular 2s 1/ / 2 state since 
there is a significant probability of the neutron in 11 Be being near the 10 Be 
core. Also, comparisons with full finite range calculations in the DWBA and 
CDCC cases, where the n-p interaction was taken to be a central Hulthcn 
potential, showed this to be a reasonable first approximation. 

The entrance and exit channel potentials are obtained from the global pa- 
rameterization of Bechetti and Greenless . We note however at this point 
that for the rigorous reproduction of experimental data, one needs to consider 
different combinations of optical potentials for the proton and deuteron chan- 
nels in the calculations in order to test the sensitivity of extracted spectroscopic 
factors to the input parameters. Despite the uncertainties on optical potentials 
used in our calculations, which in fact are important for a precise description 
of transfer cross sections, the results can however be expected to describe re- 
alistically the dynamical features of such reactions. A clear improvement for 
future work and the analysis of data would be achieved by measuring elastic 
scattering together with transfer reactions such that empirical information on 
optical potentials is obtained. 

The radial integrals are carried out from to 35 fm in steps of 0.1 fin. 
The maximum number of partial waves used is 30 for both entrance and exit 
channels. The transferred neutron wavefunction is evaluated in a Woods-Saxon 
well with shape parameters r=1.25 fm and a=0.65 fm. The real well depth is 
adjusted to reproduce the neutron separation energy. The spin-orbit force in 
the proton channel is fixed at 6 MeV. The spectroscopic factors are set to unity 
throughout the calculations. 




(25) 
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All calculations presented here are done without non-locality corrections. 
Such corrections for halo transfer are expected to be small || because they 
correct the transition amplitude in the nuclear interior, but the long tail of the 
halo wavefunction makes internal contributions less important. 



5 Results And Discussion 

From the results shown in Figs.?? and ??, it is clear that the inclusion of 
deuteron breakup is of importance for extracting information on the + and 2 + 
states. The calculated cross-section angular distributions also show significant 
differences when using the different breakup treatments. In addition, cross- 
sections of the 11 Be(p,d) 10 Be (g.s., + ) reaction calculated with different three- 
body models are larger but decrease faster at forward angles when compared 
with the calculations using DWBA which do not account for the effects arising 
from the breakup of the deuteron in the field of the nucleus. To provide an 
assessment of the three-body models (including spin-orbit interactions) used 
here we have compared in Fig. ?? the calculations (for the 2sx/a transition) 
with those obtained using the CDCC technique (excluding spin-orbit forces). 
Based on this comparison, the AD, QAD and POA are reliable tools for the 
spectroscopic study of the reaction of interest. The angular distributions of the 
transfer reactions to the 2 + state of 10 Be at small angles look the same both 
with and without breakup effects but have different absolute values. 

Our calculations have demonstrated that the inclusion of breakup effects in 
transfer reactions produces changes in the shape of the angular distributions 
and these effects increase the absolute values of the theoretical cross-sections 
in forward direction and thus lead to the smaller values of the spectroscopic 
factors extracted from the experimental data. It is expected that such effects 
will be larger when the incident energy increases and the mass of the target 
decreases f§,||. 

We therefore conclude that the standard procedure for the determination 
of spectroscopic factors as ratios of the experimental transfer cross-sections to 
those calculated within the standard DWBA is not reliable for the reactions in- 
volving weakly bound halo nuclei. The spectroscopic factors extracted from the 
experimental data would be smaller than those obtained with the conventional 
DWBA, due to including breakup effects. 

Our results also suggest that deuteron breakup stronger for 2s]/2 transfer 
than for lds/2 transfer. This may be associated with the node in the bound state 
wavefunction in the 2si/2 case. In the case of the 2si/2 transfer the shape of the 
differential cross-sections changes more strongly at small angles compared to 
the ldg/2 transfer. This should influence the ratio of the spectroscopic factors 
for + and 2 + states of the 10 Be obtained with different theoretical models. 

Fig. ?? shows the ratio of the calculated cross-sections relevant to deducing 
only relative spectroscopic information. The figure leads us to the conclusion 
that ratios of spectroscopic factors depend on possible uncertainties on abso- 
lute cross-section values, and are typically dependent upon the ingredients of 
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reaction calculations with respect to scattering angles. 

In the broad field of today's nuclear structure research , there is still very 
much to be understood when we analyze light exotic nuclei. Our intention for 
the near future is to have a better insight into the physical nature of the less 
known halo systems, using the models employed in this work. In particular, the 
question about the halo nature of 19 C and its underlying structure is one of 
the interesting current questions in driplinc physics. More information on this 
carbon isotope will allow us to further explore the characteristics of the halo 
phenomenon and to test the concepts developed for 11 Be. 
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Figure Captions 

Fig ?? Calculated differential cross-section angular distributions for the 
11 Be(p,d) 10 Be (g.s.) reaction at 35 McV using different theoretical models. 
The spectroscopic factor is 1.0 for all calculations. 

Fig. ?? As for Fig. ??, but for the n Be(p,d) 10 Bc (2+, 3.368 McV) state 
transition. 

Fig. ?? The ratio of the calculated relative cross-sections to the 2 + and + 
10 Bc final states. 



